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We find the static multi-soliton solutions of the baby Skyrme model on the two-sphere for topolog- 
ical charges 1 < i? < 14. Numerical full-field results show that the charge-one Skyrmion is spherical, 
the charge-two Skyrmion is toroidal, and Skyrmions with higher charge all have point symmetries 
which are subgroups of 0(3). We find that a rational map ansatz yields very good approximations 
to the full-field solutions. We point out a strong connection between the discrete symmetries of our 
solutions and those of corresponding solutions of the 3-D Skyrme model. 
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I. INTRODUCTION 

The Skyrme model [l| is a non-linear theory of pions 
in (3-1-1) dimensions with topological soliton solutions 
called Skyrmions. The existence of stable solutions in 
this model is a consequence of the nontrivial topology of 
the mapping M of the physical space into the field space 
at a given time, M : ^ SU{2) = 5*^, where 
the physical space 'E? is compactified to by requir- 
ing the spatial infinity to be equivalent in each direction. 
The topology which stems from this one-point compacti- 
fication allows the classification of maps into equivalence 
classes, each of which has a unique conserved quantity 
called the topological charge. 

Although Skyrmions were originally introduced to de- 
scribe baryons in three spatial dimensions [1], they have 
been shown to exist for a very wide class of geometries 
Q , and are now playing an increasing role in other areas 
of physics as well. 

Apart from the original 3D model which exhibits very 
structured solitonic solutions (for a review see ^ ) , other 
Skyrme models are known to yield solutions with intri- 
cate structures. The baby Skyrme model, first introduced 
by [3| , is a two dimensional version of the original Skyrme 
model, with as its domain. As its older brother, it is 
known to give rise, under certain settings, to structured 
multi-Skyrmion configurations [1, H, H, 0, Q • Studies of 
other Skyrme models defined on curved domains, such 
as two- and three-spheres can also be found in the lit- 
erature 0, [Hi ■ Although most of these models 
are used as a simplification or as 'toy' models of the full 
iD model, they also have physical significance on their 
own, having several potential applications in condensed- 
matter physics while often revealing useful mathematical 
features. 

In the present paper we consider a baby Skyrme model 
on the two-sphere. This type of model has been studied 
in 0, [loj , where only rotationally-symmetric configura- 



tions have been considered. We compute the full-field 
minimal energy solutions of the model up to charge 14 
and show that they exhibit complex multi-Skyrmion solu- 
tions closely related to the Skyrmion solutions of the 'iD 
model with the same topological charge. To obtain the 
minimum energy configurations, we apply two completely 
different methods. One is a full-field relaxation method, 
with which exact numerical solutions of the model are 
obtained. The other approach is a rational map approx- 
imation scheme, which as we show yields very good ap- 
proximate solutions. We discuss these methods in detail 
in section Hill 

In an exact analogy to the 31? Skyrme model, our re- 
sults show that the charge-one Skyrmion has a spherical 
energy distribution, the charge-two Skyrmion is toroidal, 
and Skyrmions with higher charge all have point sym- 
metries which are subgroups of 0(3). The symmetries 
of these solutions are the same as those of the 3D 
Skyrmions. As we shall see, this is not a coincidence. 



II. THE BABY SKYRME MODEL ON THE 
TWO-SPHERE 

The model in question is a baby Skyrme model in 
which both the domain and target are two-spheres. It 
consists of a triplet of real scalar fields (p = {4>i, 4>2,'p3) 
subject to the constraint (f> ■ (p = 1- The Lagrangian 
density is simply 
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+ - [{d^<f> ■ d^^cpf - (9^0 • dA){d^(p ■ d''<l>)\ , 

with metric ds^ = c?t'^— d^'^— sin'^ 9 d<p^, where 6 is the po- 
lar angle G [0,7r] and Lp is the azimuthal angle G [0, 27r). 
The Lagrangian of this model is invariant under rota- 
tions in both domain and the target spaces, possessing 
an 0(3)dojnain ^ 0(3)target symmetry. The first term in 
the Lagrangian is the kinetic term and the second term, 
of the fourth order in derivatives, is the 2D analogue 
of the Skyrme term [4|. While in flat two dimensional 
space a third potential term is necessary to ensure the 
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existence of stable solutions, in the present model it is 
not. Furthermore, stable solutions exist even with the 
Skyrme term missing;. This is the well known 0(3) (or 
CP^) sigma model jl3| . 

The field 4> in this model is an S*^ mapping. The 

relevant homotopy group of this model is tt2{S'^) = Z, 
which implies that each field configuration is character- 
ized by an integer topological charge B, the topological 
degree of the map 0, which in spherical coordinates is 
given by 
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where dfi — ainO AO ALp. 

Static solutions within each topological sector are ob- 
tained by minimizing the energy functional 
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Before proceeding, it is worth noting here that setting 
K = in Eq. ^ yields the energy functional of the 0(3) 
sigma model. The latter has analytic minimal energy 
solutions within every topological sector, given by 

= (sin f{e) coa{Bip), sin f{e) am{Bip), cos f{e)) , (4) 

where f{e) = cos-i(l - 2(1 -f- _(Atan6'/2)2S)-i) with A 
being some positive number [l3|. These solutions are 
not unique, as other solutions with the same energy may 
be obtained by rotating (jH) either in the target or in 
the domain spaces. The energy distributions of these 
solutions in each sector are rotationally symmetric, with 
total energy Eb = AttB. 

We have found that analytic solutions also exist for 
the energy functional ([3]) with the Skyrme term only. 
They too have the rotationally symmetric form ^ with 
/(0) = 6 and total energy Eb — AttB^. They can be 
shown to be the global minima by the following Cauchy- 
Schwartz inequality: 
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(5) 



The LHS is simply B^ and the first term in parenthesis on 
the RHS integrates to 1 . Noting that the second term in 
the RHS is the Skyrme energy (without the k^/2 factor), 
the inequality reads E > AnB^, with equality for the 
rotationally-symmetric solutions . 



III. 



STATIC SOLUTIONS 



In general, if both the kinetic and Skyrme terms are 
present, static solutions of the model cannot be obtained 



analytically. The Euler-Lagrange equations derived from 
the energy functional ^ are non-linear PDE^s, so the 
minimal energy configurations can only be obtained with 
the aid of numerical techniques. This is with the ex- 
ception of the B = I Skyrmion which has an analytic 
"hedgehog" solution 

= (sin cos sin 61 sin (p, cos 6*) , (6) 

E 

with total energy — = 1 H . 

An 2 

For Skyrmions with higher charge, we find the minimal 
energy configurations by utilizing a full-field relaxation 
method, described in more detail below. In parallel, we 
also apply the rational map approximation method, orig- 
inally developed for the 3-D Skyrme model and directly 
compare the results with the relaxation method. This 
method is also discussed below. 



A. Full-field relaxation method 

For the relaxation method, the domain S"^ is dis- 
cretized to a spherical grid - 100 grid points for 9 and 
100 points for Lp. The relaxation process begins by ini- 
tializing the field triplet to a rotationally-symmetric 
configuration 

"^initial ~ (^^^ ^ '^'^'^ ^^'^ ^ ^) ' 

where B is the topological charge of the Skyrmion in 
question. The energy of the baby Skyrmion is then min- 
imized by repeating the following steps: a point {9„n fn) 
on the grid is chosen at random, along with one of the 
three components of the field (j){9m, fn)- The chosen 
component is then shifted by a value 6^ chosen uniformly 
from the segment [— A0,A0] where = 0.1 initially. 
The field triplet is then normalized and the change in en- 
ergy is calculated. If the energy decreases, the modifica- 
tion of the field is accepted and otherwise it is discarded. 
The procedure is repeated while the value of A^ is grad- 
ually decreased throughout the procedure. This is done 
until no further decrease in energy is observed. 

One undesired feature of this minimization scheme is 
that it can get stuck at a local minimum. This prob- 
lem can be resolved by using the "simulated annealing" 
algorithm [13, [3]) which in fact has been successfully im- 
plemented before, in obtaining the minimal energy con- 
figurations of static two and three dimensional Skyrmions 
^16i] . The algorithm comprises of repeated applications of 
a Metropolis algorithm with a gradually decreasing tem- 
perature, based on the fact that when a physical system 
is slowly cooled down, reaching thermal equilibrium at 
each temperature, it will end up in its ground state. This 
algorithm, however, is much more expensive in terms of 
computer time. We therefore employ it only in part, just 
as a check on our results, which correspond to a Metropo- 
lis algorithm of zero temperature. 
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As a further verification, we set up the minimization 
scheme using different initial configurations and grids of 
different sizes (80 x 80 and 200 x 200) for several k and 
B values. This was done to make sure that the final 
configurations are independent of the discretization and 
cooling scheme. Accuracy was also verified by checking 
conservation of the topological charge B throughout the 



minimization process, yielding 
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< 10" 



B. The rational map ansatz 

Computing the minimum energy configurations using 
the full non-linear energy functional is a procedure which 
is both time-consuming and resource-hungry. To cir- 
cumvent these problems, the rational map ansatz scheme 
has been devised. First introduced in this scheme 
has been used in obtaining approximate solutions to the 
Skyrme model using rational maps between Riemann 
spheres. Although this representation is not exact, it 
drastically reduces the number of degrees of freedom in 
the problem, allowing computations to take place in a 
relatively short amount of time. The results in the case 
of 31? Skyrme model are known to be quite accurate. 

Application of the approximation, begins with express- 
ing points on the base sphere by the Riemann sphere 

Q 

coordinate z = tan— e''''. The complex- valued function 

R{z) is a rational map of degree B between Riemann 
spheres 



i?(z) 
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where p(z) and q{z) are polynomials in z, such that 
max[deg(p), deg((7)] = B, and p and q have no common 
factors. Given such a rational map, the ansatz for the 
field triplet is 



R-R 1-\R[ 
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It can be shown that rational maps of degree B corre- 
spond to field configurations with charge B [13] . Sub- 



stitution of the ansatz ([9]) into the energy functional Q 
results in the simple expression 



with 
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Minimizing the energy (jlOp only requires finding the ra- 
tional map which minimizes the functional I. As we 
discuss in the next section, the expression for X given in 
Eq. (jlip is encountered in the application of the rational 
map in the context of Skyrmions, where the proce- 
dure of minimizing I over all rational maps of the various 
degrees has been used [itI, [H, [IH . 

Here we redo the calculations, using a relaxation 
method. To obtain the rational map of degree B that 
minimizes X, we start off with a rational map of degree 
B, with the real and imaginary parts of the coefhcients of 
p(z) and q{z) assigned random values from the segment 
[—1,1]. As in the full- field relaxation method discussed 
above, solutions are obtained by relaxing the map until 
a minimum of X is reached. 



IV. RELATION TO THE 3D SKYRME MODEL 

In the 3D Skyrme model, the rational map ansatz can 
be thought of as taken in two steps. First, the radial 
coordinate is separated from the angular coordinates by 
taking the SU(2) Skyrme field U{r, 6, (p) to be of the form 



U{r,0,^)^exp{if{r) ct,{d,^) 



(12) 



where cr = (di, cr2, (T3) are Pauli matrices, f{r) is the 
radial profile function subject to the boundary conditions 
/(O) ^ n and f{oo) = 0, and cf>{9,ip) : ^ is a 
normalized vector which carries the angular dependence 
of the field. In terms of the ansatz (fT^ . the energy of 
the Skyrme field is 



E = jAnfr^Ar + J 2(/' V 1) sin^dr / ((a,</.)^ + -^(9,0)^) dl7 + J ^dr J (13) 

I 



Note that the energy functional is actually the en- 
ergy functional of our model ([3]) once the radial coor- 
dinate is integrated out. Thus, our 2D model can be 
thought of as a 3D Skyrme model with a 'frozen' radial 
coordinate. 



The essence of the rational map approximation is the 
assumption that </)(0, if) takes the rational map form ([9]), 
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which in turn leads to a simple expression for the energy 

E^AttJ (^rV' + 2B{r' + 1) sin^ / + dr , 

(14) 

where X is given in Eq. pip . As in our case, minimizing 
the energy functional requires (as a first step, followed 
by finding the profile function /(r)) minimizing X over 
all maps of degree B. 

Since the symmetries of the 3D Skyrmions are deter- 
mined solely by the angular dependence of the Skyrme 
field, it should not be too surprising that the solutions 
of the model discussed here share the symmetries of the 
corresponding solutions of the 3D Skyrme model. 



V. RESULTS 

The two approaches discussed in section IIIII have 
been applied separately to obtain the static solutions 
for charges 2 < B < lA (the charge-one solution has 
an analytic representation, as discussed in section 
This was done for several k values within the range 
0.01 < < 0.2, although solutions with different k's 
are qualitatively similar. 

As discussed in the previous section, the configura- 
tions obtained from the full-field relaxation method were 
found to have the same symmetries as corresponding 
multi- Skyrmions of the 3D model with the same charge. 
The B = 2 solution turned out to be axially symmet- 
ric, whereas higher-charge solutions were all found to 
have point symmetries which are subgroups of 0(3). For 
B = 3 and B — 12, the Skyrmions have a tetrahedral 
symmetry. The B — A and B = 13 Skyrmions have a 
cubic symmetry, and the i? = 7 is dodecahedral. The 
other Skyrmion solutions have dihedral symmetries. For 
B = 5 and B = 14 a D2d symmetry, for B = 6,9 and 
10 a D^d symmetry, for B = 8 a D^d symmetry and for 
_B = 11 a D^h symmetry. The energy distributions of the 
obtained solutions for — 0.05 are shown in FIG. [1] 

While for solutions with B < 8 the energy density (and 
also the charge density) is distributed in distinct peaks, 
for solutions with higher charge it is spread in a much 
more complicated manner. 

The total energies of the solutions (divided by AttB) 
are listed in Table [H along with the symmetries of the 
solutions (again with — 0.05). 

Application of the rational map ansatz yielded results 
with only slightly higher energies, only about 0.3% to 
3% above the full-field results. The calculated values of 
X were found to agree with those obtained in [l^ in the 
context of 3-D Skyrmions. For 9 < B < 14, the rational 
map approximation yielded slightly less symmetric solu- 
tions than the full-field ones. Considering the relatively 
small number of degrees of freedom, this method all-in- 
all yields very good approximations. The total energies 
of the solutions obtained with the rational map approxi- 
mation is also listed in Table HI 



TABLE I: Total energies (divided by 47rB) of the 
multi-solitons of the model for — 0.05. 

Charge Total energy Total energy DiflFerence Symmetry of 



B 


Full-field 


Rational maps 


in % 


the solution 


2 


1.071 


1.073 


0.177 


Toroidal 


3 


1.105 


1.113 


0.750 


Tetrahedral 


4 


1.125 


1.129 


0.359 


Cubic 


5 


1.168 


1.179 


0.958 


Did 


6 


1.194 


1.211 


1.426 


Did 


7 


1.209 


1.217 


0.649 


Icosahedral 


8 


1.250 


1.268 


1.406 


Dm 


9 


1.281 


1.304 


1.771 


Did 


10 


1.306 


1.332 


1.991 


Did 


11 


1.337 


1.366 


2.224 


Dih 


12 


1.360 


1.388 


2.072 


Tetrahedral 


13 


1.386 


1.415 


2.137 


Cubic 


14 


1.421 


1.459 


2.712 


D2 



VI. SUMMARY AND CONCLUSION 

We have studied the baby Skyrme model on the two- 
sphere, obtaining the minimal energy configurations for 
all charges up to B = 14, using both the full-field re- 
laxation method and the rational map approximation 
scheme. For each charge we have identified the sym- 
metry and measured the energy of the minimal energy 
Skyrmion. The solutions turned out to yield very struc- 
tured configurations, exhibiting the same symmetries as 
the corresponding solutions of the 3D Skyrme model. 

We have explained the reason for these similarities be- 
tween the symmetries of the two models and in the pro- 
cess we have exhibited a strong connection between them. 
The model discussed in this paper may be thought of 
as the 3D Skyrme model with a 'frozen' radial coordi- 
nate. In that sense, our computations may serve as an 
additional corroboration of results obtained for the 3D 
Skyrmions. 

In addition, we have shown that the rational map 
ansatz provides a very good approximate description to 
the true solutions, also for high topological charges. The 
energies of the solutions computed in the rational maps 
approximation are only slightly higher than the full-field 
solutions, and their symmetries, in most cases, are the 
same. This suggests that rational maps may be used to 
construct good approximations to multi- Skyrmion solu- 
tions in this model in a rather simple way. 

We believe that this work may provide a useful tool in 
the study of 3D Skyrmions, as it shares great similarities 
with the 3D model, especially in terms of multi-Skyrmion 
symmetries. The fact that the model discussed here is 
two-dimensional makes it simpler to study and perform 
computations with, when compared with the 3D Skyrme 
model. 
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(m) S = 14 

FIG. 1: The energy distributions of the multi-Skyrmion solutions for charges 2 < B < 14 (k^ = 0.05). 
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